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ABSTRACT 

The purpose of this paper is to present a method to self-calibrate the spectral 
energy distribution (SED) of objects in a survey based on the fitting of an 
SED library to the observed multi-color photometry. We adopt for illustrative 
purposes the Vilnius (Strizyz and Sviderskiene 1972) and Gunn & Stryker 
(1983) SED libraries. The self-calibration technique can improve the quality of 
observations which are not taken under perfectly photometric conditions. The 
more passbands used for the photometry, the better the results. This technique 
has been applied to the BATC 15-passband CCD survey. 

Subject headings: survey, techniques: photometric, stars: fundamental 
parameters 
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1. Introduction 

Multi-color photometry can provide accurate SED information of low spectral 
resolution for a large sample of objects and is a very powerful tool to deal with many 
important astrophysical problems. It can be used to measure the photometric red-shift of 
quasars and galaxies, and to select quasars and other interesting objects based on their 
characterstic SEDs. One can also use SEDs for stellar population synthesis, an important 
tool for studying the structure and evolution of galaxies. The more passbands in which one 
observes, the better one can determine the SED. 

The Beijing- Arizona- Taipei-Connecticut (BATC) multicolor photometric survey is 
designed to obtain as much information on the SED of celestial sources possible. It 
combines the Beijing Astronomical Observatory's 60/90 cm f/3 Schmidt telescope with 
a Ford Aerospace CCD of 2048x2048 pixels (recently replaced by a thinned Loral CCD, 
thanks to the Steward Observatory) and 15 intermediate-band filters ranging from 320 
nm to 1000 nm to obtain the SEDs for all objects down to B — 21 mag in a field of 
one square degree. A key problem in the survey is calibrate accurately the SED of our 
objects. The standard procedure of SED calibration is as follows: on photometric nights, 
we observe through the 15 BATC filters both Oke & Gunn (1983) spectroscopic standard 
stars (HD19445, HD84937, BD+26°2606 and BD+17°4708), and the target fields. We 
convolve the known fluxes of the standard stars with the transmission curves of the BATC 
filters to determine magnitudes in BATC photometric system (Fan et al. 1996, Zheng al. 
1999). In practice, however, this method is not very efficient, for two reasons. First, the 
calibrations are not always taken on perfectly photometric nights. This results in systematic 
troughs or bumps in the SED for all objects in a field. Second, in order to obtain accurate 
SED calibration for 15 passbands, we need many photometric nights, which could take a 
very long time at our present site. In order to solve these problems, we are developing a 
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technique we call "SED self-calibration" . This is a statistical method, based on applying a 
stellar SED library to a large sample of observed stellar objects. In section 2, we describe 
our method of SED self-calibration. In section 3 we present some results from tests of the 
method. We discussions the method and give our conclusions in section 4. 



2. Method 

The BATC survey covers regions at high galactic latitudes. In a typical field, several 
thousand objects are detected per square degree. Among them are several hundred 
unsaturated, bright stars with high signal-to-noise ratios. These "good" stars of highly 
precise photometry form a big sample which is the observational basis for the SED 
auto-calibration method. In order to make the description concise, we will use, in the 
following presentation, the term SED obs to express the SED of an object which is calibrated 
by observations, and the term SED matc h to express the SED of the same object which is 
the closest match found in the SED library. We assume that most of these "good" stars 
are normal, nearby stars; therefore, most of these "good" stars should have very close 
matches in the SED libraries. If the field has been well calibrated by observations, the 
RMS of the differences between the SED obs and the SED match for most of the "good" stars 
should be roughly the size of the photometric precision. If the RMS of the residuals is 
considerably larger than the photometric precision, then either some of the "good" stars 
might not actually be normal stars, or observations in some passbands might not have been 
made during truly photometric nights. In the former case, we can reject those stars with 
abnormal residuals from the data sample and repeat the process. In the later case, we can 
shift the zero points for some passbands which show systematic deviations of SED obs from 
SED matc h- We iterate these processes until the RMS of the residual reaches the level of 
the photometric precision, which is about 0.01 mag for "good" stars. This self-calibration 
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method is based on the following assumptions: 

1. Most bright stellar objects in the field are normal stars whose SEDs can be found in 
the SED library. 

2. The interstellar extinction of these bright stars is small and the differences in 
extinction between stars in the same CCD field is negligible. 

3. The SED library is reliable and covers all the required types of spectral and 
luminosity classes to fit the SED obs . 

We now describe briefly details of the observed and model spectra energy distributions. 

2.1. SED bs 

First, we must obtain SED obs for all the objects in the target field. Once a target 
field has been observed in several passbands, we can obtain the instrumental colors of the 
stars in the field by using standard photometry packages, such as DAOPHOT. The zero 
point for each passband is obtained through observations of Oke & Gunn standards during 
photometric nights. As long as the zero points are well-determined, the instrumental 
magnitudes can be transferred into the calibrated magnitudes of the BATC system. The 
instrumental color index can be defined as in Table 1, where Co] . is the instrument color 
index of j'th band minus s'th band of the i'th object. Here we use the s'th band as the 
reference band 

To calculate the flux-calibrated color index (7j , the color index zero point correction 
Ccj tS is added to the instrumental color index Cj tS . 



CI = Col + Co,, 
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Table 1: Table of the instrumental color index 

2.2. SED libraries 

The stellar SED library used for this method is a hybrid one, including the theoretical 
SEDs of Kurucz (1992, 1993), and the observational SEDs of Gunn and Stryker, and of 
Straizyz and Sviderskiene (henceforth referred to as "Vilnius"). We transfer spectral fluxes 
into the BATC photometric system using the following equation: 

m % = -2.5log(f v )i - 48.6 

where (f u )i is the monochromatic energy at the central wavelength of i th filter, in unit 
of ergcm~ 2 s~ 1 Hz~ 1 (cf. Fan 1995; Fan et al. 1996). The Vilnius library has 49 spectra 
covering spectral type from O to M6 and luminosity from main sequence to giant. The 
library of Gunn and Stryker contains 74 spectra with the same coverage of the spectral 
type and luminosity. The SED libraries given by Kurucz cover temperatures from 3500K to 
50000K with a range of stellar surface gravities and metallicities. Using these SEDs we can 
build a table of model color indices: 

If the SEDobs is well determined by the observations, it can be matched to one of 
the model SEDs with a residual at the level of photometric precision. Mathematically, we 
search for a minimum of a: 
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Table 2: Table of the model SED, where Cm!- is the color index between the BATC i'th 
and the s'th filter bands of the SED library for the i'th spectral type and luminosity 



a — J2 m i n 



^(Co^ + Cc^-Cm^) 2 ,fcel, 



N 



Here, n is the number of the "good" stars, m is the number of the color index, and 
N is the total number of entries in the model SED library. Co* s , Ccj :S and Cm^ s are the 
instrumental color index, the color index zero-point offset, and the color index of model 
SED, respectively. We look simultaneously for the nearest match to each star in the library, 
and a value for Cc which leads to a minimum in a. In our algorithm, we call subroutines of 
the MINUIT package (James, 1994) from the CERNLIB software. The process converges on 
the The correct values of the color index zero-point corrections Ccj yS . During the iteration 
process, a different weight is given to each star according to its instrumental magnitude. 
The program can subsequently reject those stars of having abnormally high differences 
between SED obs and SED match . At the start of the iteration, we must provide initial values 
for the color index corrections. The effect of different initial values on the final converged 
result is less than 0.01 mag. We take the mean instrumental color index and the mean 
model color index as the initial corrections: 
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Because most nearby stars are spectral type F, G and K, we employ only these types 
of model SEDs to estimate the initial color correction constants. 

3. Testing 

3.1. Empirical comparison of the two methods 

We can test the method as follows: if we observe one of the Oke & Gunn standard star 
fields, we acquire the instrumental SED for the standard star. Since we know the real SED 
of the standard star, we can determine exactly the SED corrections for each passband. In 
other words, for this particular field, the zero-point color corrections can be derived directly. 
If we use our method on the same data set, we can compare its values for the zero-point 
corrections to the correct ones. 

We observed the field of HD84937 through 13 filters on Jan. 22, 1998. The transparency 
was good, but it was not photometric. Two exposures were taken for each filter: a short 
one of a few seconds to avoid saturating the bright star HD84937, and a long one of 300 
seconds. The short and the long exposures for each filter were taken in quick succession 
in order to guarantee that both shared the same weather conditions. The short exposure 
was used to determine the SED corrections for the field via HD84937 directly, and the long 
exposure was used to determine the SED corrections via our method. In the following, we 
use this data set to do several tests of our method. 

As mentioned above, we obtained the instrumental SED of standard star HD83927 
(M AG stdJnstr) from the short exposure images. Using the known SED of HD83927 in 
BATC system (M AG std-BATC) , we calculated the differences in each passband v. 

dMAGstdi = MAGstd_BATCi - MAGstdAnstn 
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If the SED self-calibration method is successful, it should yield corrections for the SED 
of the long exposure images Cq which differ from dM AG _stdi only by a constant K, due 
to the difference in exposure time: 

K = Cd- dMAGstdi 

The results are listed in Table 3. It shows that the constant K is indeed the same in 
each passband, except for the BATC2 filter. The RMS of the variation of K in color is of 
the level of 0.004. We give two reasons for the large difference in the near-UV band BATC2: 
first, our thick CCD has very low quantum efficiency in the blue, causing low signal-to-noise 
values in all stars. Second, the SED templet in this region of wavelengths is available only 
for some stellar types. 

3.2. Cross checking between two spectral libraries 

We can test our method in a second way. Suppose the SED libraries, either Gunn 
& Stryker or Vilnius, represent well the SED for most types of stars. Having computed 
synthetically the BATC magnitudes in all passbands for the stars in each library, we can 
apply our method to calculate the corrections between the two sets of spectral models. If: 

1) both libraries provide accurate stellar SEDs, and 

2) each library covers an equal range of stellar types, i.e. any SED in one library has a 
corresponding SED in the other library, and 

3) the method developed here is correct and effective, 

then after the iteration process, the final SED corrections should be close to zero. 
In the following test, we use Gunn & Stryker library as SED match and the Vilnius 
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Table 3: Results of a test of the SED self-calibration method, using long and short exposures 
for the field of standard star HD84937 

library as the SED ohs catalog. The results are shown in Table 4. Column 1 is the filter 
name, column 2 its central wavelength, Column 3 the correction for each band. The BATC9 
passband (6660A) is used as the reference band, and kept fixed during the iterative process. 
The final values for the zero point corrections are indeed very small. The mean deviation of 
zero point corrections is ±0.015 mag; however, the BATC1 band again shows much larger 
devaitions than the other bands. 

This result shows that the two SED libraries in general can be matched well each other, 
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but two points call for for further discussion. 

1) There is a systematic deviation for both SED libraries in the blue versus the red: 
below 455nm, the Vilnius SEDs are flatter than those of Gunn & Stryker; above 455nm, 
the Vilnius SEDs are more depressed those of Gunn & Stryker. 

2) The large deviation in UV band indicates that either one or both the SEDs does not 
well represent the stellar SEDs in this region. 
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Table 4: Results of cross checking between Vilnius and Gunn & Stryker model spectra 
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Fig. 1. — Spectral type of Vilnius vs spectral types obtained from Gunn & Stryker model 

3.3. Stellar classification 

The process of matching SED obs to SED match produces as a byproduct a spectral 
classification for each star; or, if a library of theoretical models is used, the stellar 
atmospheric parameters of each star. It provides us with an indirect way to test the method 
presented in this paper. Figure 1 shows a good correlation between the original spectral 
class of stars in the Vilnius catalog, and the spectral class determined for those stars via 
our method, using the Gunn & Stryker library for SED match . 
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4. Discussions and conclusions 

1. Our method is a statistical one, so it can only be effectively applied to a large 
stellar sample. The BATC survey satisfies this condition. Each BATC field covers about 
one square degree. Typically, there are more than 4000 objects detected in each image. 
We are guaranteed several hundred bright and unsaturated stars with reliable instrumental 
magnitudes to use as a "good" star sample. 

2. The basis of our method is fitting the observed stellar SED to a library of stellar 
SEDs. Our assumption that most of our "good" stars are normal stars appears firm, and 
their SEDs can be found from the SED library. Most of abnormal stars are rejected during 
of the iteration process. 

3. The SED library is a key for the method. There is a great demand in the 
astronomical community for reliable SED libraries which cover a wide range in wavelength 
and stellar type, The theoretical libraries (e.g. Kurucz 1993) are limited by our knowledge 
of stellar physics. The theoretical SEDs for late type, low-temperature stars is not as 
reliable as those for hot stars. There are many observational stellar libraries (Gunn & 
Stryker, Vilnius, etc), but differences among them exist. It is not easy to judge which 
one is the best. After comparing various stellar SED libraries we have collected from the 
literature, we are at present using mostly the Gunn & Stryker library in our calculations. 
Further work on making good libraries of SEDs is necessary. 

4. In principle, interstellar extinction must be considered in the final result. Our 
method of SED correction takes into account not only terrestrial atmospheric extinction, but 
also interstellar reddening. If standard stars or spectral models suffer from systematically 
more or less (or different) extinction than a survey's target stars, our method will yield 
systematically incorrect zero-point corrections to the photometry of target stars. This 

is because the method uses a fitting process to fit the observational SED affected by 
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the reddening to the theoretical SED unaffected by the reddening. In the BATC survey, 
the fields are located at high galactic latitudes, where the interstellar extinction is much 
reduced. We therefore use only bright stars, which are nearby and little affected by 
interstellar extinction, as standards for the survey fields. 

5. We have made several tests to see if the method works. 

(1) We have tried different filter bands as the reference band to see if the change in the 
reference band affects the results. Our results show that the difference in the final constant 
corrections is less than 0.01 mag. This means that the choice of reference band in color 
index is not important. Normally, we select the band of deepest exposure as the reference 
band, since it has highest signal to noise ratio. This also allows the measurement of color 
index for as many stars as possible. 

(2) We wanted to know if the number of the filters used affects the results, and 
what is the minimum number of colors required to make the SED self-calibration work 
effectively. In order to do the test, we reduced the number of filters by taking off some of 
the color indexes. We found that the method still yields reasonable results. Obviously, the 
results improve when one uses a larger number of the filter bands, and a wider range of 
wavelengths. 

(3) . We wanted to know if the final results depend on the characteristics of the 
"good" stars, which we used for the SED self-calibration. For this purpose, we used various 
randomly selected subsets of the "good" stars in each image. We obtained very similar 
color correction constants. Furthermore, we divided the "good" stars sample into several 
sub-groups according to their apparent magnitude. The results show that the difference 
among different groups is about 0.03 mag. One possible reason is the low signal-to-noise 
ratio of faint stars. 
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(4) Does the choice of the initial value of the iteration affect the results? Our tests 
show that different initial values normally only affect the convergence time of iteration. The 
difference on the final result is less than 0.01 mag. In some cases, there are several minima 
and the process may not iterate to the right one; this is an open issue. The best way to 
solve this problem, from our experiences, is that if one or a few color indexes are very well 
determined by observations taken during photometric nights, then we keep these indexes 
fixed. This causes the iteration always to converge to the right minimum. 

After many tests and applications to real data, we conclude that, though there are still 
some problems requiring further development (such as creating larger libraries of stellar 
SEDs), the method presented here can work well: the accuracy of the SED calibration is 
comparable to the precision of the CCD photometry. A by-product of our method is the 
automatic classification of the stellar type or the determination of the stellar parameters, 
which is very useful for the studies of galactic structure via large field multi-color CCD 
photometry survey. 
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